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Abstract 

We present results obtained from a study of the structure of hadronic events 
recorded by the L3 detector at various centre-of-mass energies. The distributions of 
event shape variables and the energy dependence of their mean values are measured 
from 30 GeV to 189 GeV and compared with various QCD models. The energy 
dependence of the moments of event shape variables is used to test a power law 
ansatz for the non-perturbative component. We obtain a universal value of the non- 
perturbative parameter a® = 0.537 ± 0.073. From a comparison with resummed 
O(ol) QCD calculations, we determine the strong coupling constant at each of the 
selected energies. The measurements demonstrate the running of a s as expected in 
QCD with a value of a s (m z ) = 0.1215 ± 0.0012 (exp) ± 0.0061 (th). 
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1 Introduction 



LEP operated at centre-of-mass energies around 91.2 GeV from 1989 to 1995 and then moved 
up to six different centre-of-mass energies between 130 GeV and 189 GeV in the following three 
years. Thus a study of the process e + e~ — > hadrons at LEP offers a unique environment to 
test the predictions of the theory of the strong interaction (QCD) over a wide energy range. 
The energy range has been extended by using hadronic events from Z decays with isolated high 
energy photons in order to probe the structure of hadronic events at reduced centre-of-mass 
energies down to 30 GeV PUI2]- The high energy photons are radiated early in the process 
through initial state radiation (ISR) or through quark bremsstrahlung whereas the hadronic 
shower develops over a longer time scale. 

We report here measurements of event shape distributions and their moments using the data 
collected with the L3 detector |3j. We update the published results at \/s = 161, 172 and 183 
GeV jlJE] with an improved selection method for hadronic events and present new results at y/s 
= 130, 136 and 189 GeV. The measured distributions are compared with predictions from event 
generators based on an improved leading log approximation (Parton Shower models including 
QCD coherence effects). Three such Monte Carlo programs (Jetset PS jO], Herwig [Zj and 
Ariadne |Hj) have been used for these comparisons. We also compare our measurements with 
predictions from QCD models with no coherence effects (Cojets 9\). These Monte Carlo 
programs use different approaches to describe both the perturbative parton shower evolution 
and non-perturbative hadronisation processes. They have been tuned to reproduce the global 
event shape distributions and the charged particle multiplicity distribution measured at 91.2 
GeV P3j. 

The moments of event shape variables are measured between 30 GeV and 189 GeV. The 
perturbative and non-perturbative QCD contributions are obtained from a fit using the power 
correction formula [11J. This approach was first applied by the DELPHI collaboration [T2*] . 

The strong coupling constant a s is also determined at each of these centre-of-mass energies 
by comparing the measured event shape distributions with predictions of second order QCD 
calculations containing resummed leading and next-to-leading order terms [T3j . 

Section El describes the selection of hadronic events. Measurements of event shape variables 
and estimation of systematic errors are described in section El Section 0] presents a comparison 
of the data with predictions from various QCD models, a study of the power correction ansatz 
and a determination of a s from event shape distributions. The results are summarised in 
section El 

2 Event Selection 

The selection of e + e~ — > hadrons events is based on the energy measured in the electromagnetic 
and hadron calorimeters. We use energy clusters in the calorimeters with a minimum energy of 
100 MeV. We measure the total visible energy (E vis ) and the energy imbalances parallel (E\\) 
and perpendicular (E±) to the beam direction. Backgrounds are different for hadronic Z decays, 
hadronic events at reduced centre-of-mass energies and at high energies. This is reflected in 
the different selection cuts used for these three types of data sets. 

We use Monte Carlo events to estimate the efficiency of the selection criteria and purity 
of the data sample. Monte Carlo events for the process e + e~ — ► qq(7) have been generated 
by the parton shower programs Jetset and Pythia ^3] and passed through the L3 detec- 
tor simulation [TE] . The background events are simulated with appropriate event generators: 
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Pythia and Phojet jUj for two-photon events, Koralz [THj for the t + t~{j) final state, 
Bhagene [in] and Bhwide (201 f° r Bhabha events, Koralw [21] for W-pair production and 
Pythia for Z-pair production. 

Details of event selection at y/s « mz and at reduced centre-of-mass energies have been 
described earlier [HI2]- At y/s ~ m z , we have used only a small subset of the complete data 
sample (8.3 pb _1 out of 142.4 pb _1 of integrated luminosity) which still provides an experimental 
error three times smaller than theoretical uncertainties. 

Data at y/s = 130 and 136 GeV were collected in two separate runs during 1995 [I] and 1997. 
The main background at these energies comes from ISR resulting in a mass of the hadronic 
system close to m z . This background is reduced by applying a cut in the two dimensional plane 
of | E\\ | / E vis and E vis /y/s. In the current analysis, data sets from the two years have been 
combined and the cuts are optimised to get the best efficiency times purity. 

For the data at y/s > 161 GeV, additional backgrounds arise from W-pair and Z-pair 
production. A substantial fraction (~ 80%) of these events can be removed by a specific 
selection jS] based on: 

• forcing the event to a 4-jet topology using the Durham algorithm [2*2"] . 

• performing a kinematic fit imposing the constraints of energy-momentum conservation, 

• making cuts on energies of the most and the least energetic jets and on yf 4 , where yf 4 is 
the jet resolution parameter for which the event is classified as a three-jet rather than a 
four-jet event. 

These cuts have also been optimised at each energy point. For centre-of-mass energies at or 
above 130 GeV, hadronic events with ISR photon energy larger than 0.18a/s are considered as 
background. 

The integrated luminosity, selection efficiency, purity and number of selected events for each 
of the energy points are summarised in Table ^ 

3 Measurement of Event Shape Variables 

We measure five global event shape variables for which improved analytical QCD calculations 
[T4] are available. These are thrust (T), scaled heavy jet mass (p), total (B T ) and wide (B w ) 
jet broadening variables and the C-parameter. 

For Monte Carlo events, the global event shape variables are calculated before (particle 
level) and after (detector level) detector simulation. The calculation before detector simulation 
takes into account all stable charged and neutral particles. The measured distributions at 
detector level differ from the ones at particle level because of detector effects, limited acceptance 
and resolution. After subtracting the background obtained from simulations, the measured 
distributions for all energies except y/s ~ mz are corrected for detector effects, acceptance 
and resolution on a bin-by-bin basis by comparing the detector level results with the particle 
level results. The level of migration is kept at a negligible level with a bin size larger than the 
experimental resolution. At y/s « mz, the detector effects are unfolded for these event shape 
variables using a regularised unfolding method [211] • We also correct the data for initial and 
final state photon radiation bin-by-bin using Monte Carlo distributions at particle level with 
and without radiation. 
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The systematic uncertainties in the distributions of event shape variables arise mainly due 
to uncertainties in the estimation of detector correction and background estimation. The un- 
certainty in the detector correction has been estimated by several independent checks: 

• The definition of reconstructed objects used to calculate the observables has been changed. 
Instead of using only calorimetric clusters, the analysis has been repeated with objects 
obtained from a non-linear combination of energies of charged tracks and calorimetric 
clusters. At y/s ~ m z , we use a track based selection and the event shape variables are 
constructed from the tracks. 

• The effect of different particle densities in correcting the measured distribution has been 
estimated by changing the signal Monte Carlo program (Herwig instead of Jetset). 

• The acceptance has been reduced by restricting the events to the central part of the 
detector (| cos(#t)| < 0.7, where # T is the polar angle of the thrust axis relative to the 
beam direction) where the energy resolution is better. 

The uncertainty on the background composition of the selected event sample has been esti- 
mated differently for the three types of data sets. At « mz, the background contamination 
is negligible and the uncertainty due to that has been neglected. For data samples at reduced 
centre-of-mass energies, the systematic errors arising from background subtraction have been 
estimated [2^ by: 

• varying, by one standard deviation, the background scale factor which takes into account 
the lack of isolated it and 77 production in the Monte Carlo sample, 

• varying the cuts on neural network probability, jet and local isolation angles, and energy 
in the local isolation cone. 

At high energies, the uncertainty is determined by repeating the analysis with: 

• an alternative criterion to reject the hard initial state photon events based on a cut on 
the kinematically reconstructed effective centre-of-mass energy, 

• a variation of the estimated two-photon interaction background by ± 30% and by changing 
the background Monte Carlo program (Phojet instead of Pythia), and 

• a variation of the W + W~ background estimate by changing the W-pair rejection criteria. 

The systematic uncertainties obtained from different sources are combined in quadrature. 
At high energies, uncertainties due to ISR and W + W~ backgrounds are the most important 
ones. They are roughly equal and are 2-3 times larger than the uncertainties due to the detector 
correction. 

Apart from the data set at y/s ~ mz, statistical fluctuations are not negligible in the 
estimation of systematic effects. The statistical component of the systematic uncertainty is 
determined by splitting the overall Monte Carlo sample into luminosity weighted sub-samples 
and treating each of these sub-samples as data. The spread in the mean position gives an 
estimate of the statistical component and is taken out from the original estimate in quadrature. 
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4 Results 



4.1 Comparison with QCD models 

Figure [T] shows the corrected distributions for thrust, scaled heavy jet mass, total and wide 
jet broadening and the C-parameter obtained at \/s = 189 GeV. The data are compared 
with predictions from QCD models Jetset PS, Herwig and Ariadne at particle level. The 
agreement is satisfactory. 

An important test of QCD models is a comparison of the energy evolution of the event shape 
variables. The energy dependence of the mean event shape variables arises mainly from two 
sources: the logarithmic energy scale dependence of a s and the power law behaviour of non- 
perturbative effects. The first moments of the five event shape variables are shown in Figure El 
and TableEJ Also shown are the energy dependences of these quantities as predicted by Jetset 
PS, Herwig, Ariadne, Cojets and Jetset ME {0(al) matrix element implementation). 
All the models with the possible exception of Jetset ME give a good description of the data. 



4.2 Power Law Correction Analysis 

The energy dependence of moments of the event shape variables has been described [IT] as a 
sum of the perturbative contributions and a power law dependence due to non-perturbative 
contributions. The first moment of an event shape variable / is written as 

(/) = (/pert) + (/pow) , (1) 

where the perturbative contribution (/ per t) has been determined to 0(a?) [21] • The power 
correction term ^JJ , for 1 — T, p, and C, is given by 

(/pow) = c f V, (2) 

where the factor Cf depends on the shape variable / and V is supposed to have a universal 
form: 



7l Z y/S 



111 \ flj Pq J 



(3) 



for a renormalisation scale fixed at \/s. The parameter «o is related to the value of a s in the 
non-perturbative region below an infrared matching scale /// (= 2 GeV); /3q is (lliV c — 2vV/)/3, 
where N c is the number of colours and Nf is the number of active flavours. K = (67/18 — 
7t 2 /6)Ca — 5Nf/9 and Cp, Ca are the usual colour factors. The Milan factor Ai is 1.49 for Nf 
= 3. For the jet broadening variables, the power correction term takes the form 

(/pow) = c f FV, (4) 

where 

( * 3 fa 



\2y/aC F a C Mw + 4 Qa C F 



0.6137 + 0(y/a^\ (5) 



and a takes a value 1 for Bt and 2 for Bw and «cmw is related to a s [TT] . 

We have carried out fits to the first moments of the five event shape variables separately with 
<^s( m z) and cto as free parameters. The diagonal terms of the covariance matrix between the 
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different energy points are constructed by summing in quadrature the systematic uncertainty 
and the statistical error. The off-diagonal terms are obtained from the common systematic 
errors. The results of the fits are summarised in Table El and shown in Figure El 

The five values of obtained from the event shape variables agree within errors, supporting 
the predicted universality of the power law behaviour. The theoretical predictions for event 
shape variables, being incomplete, give different estimates of «o and a B . Since the measurements 
are fully correlated, the best estimates of the overall values are obtained by taking an unweighted 
average: 

a = 0.537 ± 0.070 ± 0.021 , (6) 
a s (m z ) = 0.1110 ± 0.0045 ± 0.0034. (7) 

The first error on each measurement is experimental and is obtained from the average of the 
five errors on a and a s . To estimate theoretical uncertainties we vary the renormalisation 
scale between 0.5y/s and 2.0 y/s and cto and a s (mz) vary on average by ±0.021 and ±0.0033 
respectively. A variation of /xj in the range from 1 to 3 GeV gives an additional uncertainty 
on a s (mz) of ± 0.0010. These two estimates of theoretical uncertainties are combined in 
quadrature and quoted as the second error. 

We have also measured the second moments of these shape variables which are summarised 
in Table El The energy dependence of these moments has been analysed in terms of power law 
corrections. For variables 1 — T, p and C, the following result is expected to hold |25j : 

(f) = (/ P 2 crt ) + 2(/ pcrt ) C/ P + O (A . (8) 

This assumes that the non-perturbative correction to the distributions causes only a shift. 
For jet broadenings the power corrections are more complicated. The O(^) term has been 
parametrised as A 2 /s and is expected to be small for 1 — T,p and C. Fits have been performed 
to the second moments where ao and a s have been fixed to the values obtained from the 
corresponding fits to the first moments. Figure |U shows the second moments compared to these 
fits. The contributions of the 0(-) term are non- negligible for 1 — T and C, in contradiction 
with the expectation. The five values of A 2 , as obtained from the fits, are summarised in 
Table El 

4.3 ot s from Event Shape Distributions 

In order to derive a s from event shape variables at each energy point we fit the measured 
distributions to theoretical calculations based on O(ol) perturbative QCD with resummed 
leading and next-to-leading order terms. These calculations are performed at parton level 
and do not include heavy quark mass effects. To compare the analytical calculations with 
the experimental distributions, the effects of hadronisation and decays have been corrected for 
using Monte Carlo programs. 

The fit ranges used take into account the limited statistics at high energy as well as the 
reliability of the resummation calculation and are given in Table El In this analysis, we deter- 
mine ot s at y/s = 130, 136 and 189 GeV for the first time. We also include the measurements 
done at y/s = 161, 172 and 183 GeV since the experimental systematic uncertainties are con- 
siderably reduced by using an improved selection method and by subtracting the statistical 
component of the systematic uncertainties. All the measurements are summarised in Table El 
These measurements supersede those published previously jS]. 
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The experimental errors include the statistical errors and the experimental systematic uncer- 
tainties. The theoretical error is obtained from estimates [Sj of the hadronisation uncertainty 
and of the errors coming from the uncalculated higher orders in the QCD predictions. The 
estimate of the theoretical error does not always reflect the true size of uncalculated higher 
order terms. An independent estimate is obtained from a comparison of a s measurements 
from many event shape variables which are affected differently by higher order corrections and 
hadronisation effects. To obtain a combined value for the strong coupling constant we take the 
unweighted average of the five a s values. We estimate the overall theoretical error from the 
simple average of the five theoretical errors or from half of the maximum spread in the five a s 
values. Both estimates yield similar results. The combined results are summarised in Table H3 
The earlier measurements at y/s = rrtz and at reduced centre-of-mass energies determined a s 
from four event shape variables only: T, p, Bt and B\y- For comparison we also provide in 
Table |B1 the mean from these four measurements. 

We compare the energy dependence of the measured a s values with the prediction from 
QCD in Figure EK- The theoretical errors are strongly correlated between these measurements. 
The error appropriate to a measurement of the energy dependence of a s can then be considered 
to be experimental. The experimental systematic errors on a s are dominated by the back- 
ground uncertainties. These are similar for all the individual low energy or high energy data 
points but differ between the low energy, Z peak and high energy data sets. The experimental 
systematic errors are then different and uncorrelated between the three data sets, but are taken 
as fully correlated between individual low energy or high energy measurements. The thirteen 
measurements in Figure EH are shown with experimental errors only, together with a fit to the 
QCD evolution equation j2Sj with a s (mz) as a free parameter. The fit gives a x 2 of 13.5 for 12 
degrees of freedom corresponding to a confidence level of 0.34 with a fitted value of a s : 

a s (m z ) = 0.1215 ± 0.0012 ± 0.0061 . (9) 

The first error is experimental and the second error is theoretical. On the other hand, a fit 
with constant a s gives a \ 2 °f 65.1 for 12 degrees of freedom. The value of a s (mz) thus 
obtained is in agreement with the value obtained in the power law ansatz analysis considering 
the experimental and the theoretical uncertainties. 

Figure Eb summarises the a s values determined by L3 from the r lifetime measurement |27j . 
Z lineshape j2H] an d event shape distributions at various energies, together with the QCD 
prediction obtained from a fit to the event shape measurements only. These measurements 
support the energy evolution of the strong coupling constant predicted by QCD. 

The slope in the energy evolution of a s depends on the number of active flavours. We have 
performed a fit with Nf as a free parameter along with a s and obtain the number of active 
flavours: 

N f = 5.0 ± 1.3 ± 2.0 , (10) 

where the first error is experimental and the second is due to theoretical uncertainties. The 
errors have been estimated by using the covariance matrix determined from experimental and 
overall errors on a s in the fit. This result agrees with the expectation Nf = 5. 

5 Summary 

We have measured distributions of event shape variables in hadronic events from e + e~ annihi- 
lation at centre-of-mass energies from 30 GeV to 189 GeV. These distributions as well as the 
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energy dependence of their first moments are well described by parton shower models. 

The energy dependence of the first two moments has been compared with second order 
perturbative QCD with power law corrections for the non-perturbative effects. The fits of the 
five event shape variables agree with a universal power law behaviour giving a — 0.537 ± 0.070 
(exp) ± 0.021 (th). We find a non-negligible contribution from an O(^) term in describing the 
second moments of 1 — T, B T and C. 

The event shape distributions are compared to second order QCD calculations together 
with resummed leading and next-to-leading log terms. The data are well described by these 
calculations at all energies. The measurements demonstrate the running of a s as expected 
in QCD with a value of a s (m z ) = 0.1215 ± 0.0012 (exp) ± 0.0061 (th). From the energy 
dependence of a s , we determine the number of active flavours to be Nf = 5.0 ± 1.3 (exp) 
± 2.0 (th). 
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Table 1: Summary of integrated luminosity, selection efficiency sample purity and number of 
selected hadronic events at the different energies used in this analysis. The energies below 
y/s = 91 GeV are obtained from the full data sample at the Z peak, by selecting events with 
an isolated high energy photon. 
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(GeV) 


First moments of 


1 - T 


P 


Bt 


Bw 


C 


30-50 


.0971 ± .0030 ± .0034 


.0747 ± .0023 ± .0023 


.1399 ± .0027 ± .0016 


.0896 ± .0021 ± .0018 


.3667 ± .0084 ± .0073 


50-60 


.0811 ± .0027 ± .0029 


.0632 ± .0021 ± .0023 


.1223 ± .0025 ± .0054 


.0800 ± .0020 ± .0034 


.3091 ± .0080 ± .0131 


60-70 


.0796 ± .0021 ± .0051 


.0603 ± .0015 ± .0047 


.1213 ± .0019 ± .0079 


.0806 ± .0014 ± .0060 


.3049 ± .0059 ± .0232 


70-80 


.0731 ± .0015 ± .0045 


.0560 ± .0011 ± .0027 


.1157 ± .0015 ± .0048 


.0758 ± .0011 ± .0046 


.2851 ± .0044 ± .0177 


80-84 


.0700 ± .0018 ± .0046 


.0546 ± .0015 ± .0035 


.1116 ± .0017 ± .0057 


.0756 ± .0014 ± .0051 


.2759 ± .0055 ± .0191 


84-86 


.0691 ± .0022 ± .0088 


.0544 ± .0017 ± .0085 


.1102 ± .0021 ± .0086 


.0749 ± .0017 ± .0092 


.2722 ± .0068 ± .0289 


91.2 


.0636 ± .0003 ± .0013 


.0539 ± .0002 ± .0013 


.1102 ± .0002 ± .0011 


.0738 ± .0001 ± .0008 


.2599 ± .0004 ± .0054 


130 


.0556 ± .0022 ± .0014 


.0452 ± .0018 ± .0007 


.0976 ± .0023 ± .0008 


.0681 ± .0019 ± .0007 


.2277 ± .0072 ± .0052 


136 


.0614 ± .0029 ± .0011 


.0467 ± .0022 ± .0004 


.0999 ± .0029 ± .0011 


.0699 ± .0024 ± .0006 


.2357 ± .0089 ± .0038 


161 


.0513 ± .0030 ± .0008 


.0421 ± .0025 ± .0007 


.0923 ± .0032 ± .0018 


.0666 ± .0027 ± .0010 


.2052 ± .0098 ± .0028 


172 


.0542 ± .0037 ± .0022 


.0440 ± .0028 ± .0018 


.0950 ± .0046 ± .0031 


.0664 ± .0031 ± .0023 


.2281 ± .0159 ± .0133 


183 


.0539 ± .0020 ± .0011 


.0424 ± .0014 ± .0004 


.0918 ± .0020 ± .0015 


.0654 ± .0015 ± .0010 


.2157 ± .0063 ± .0073 


189 


.0548 ± .0013 ± .0013 


.0442 ± .0009 ± .0009 


.0918 ± .0013 ± .0018 


.0669 ± .0009 ± .0010 


.2160 ± .0040 ± .0041 



Table 2: First moments of the five event shape variables at different energy points. The two errors are respectively statistical and 
systematic. 



Observable 


«o 


« S K) 


xVd.o.f. 


A 2 (GeV 2 ) 


1 — T 


0.633 ± 0.097 


0.1104 ± 0.0065 


11.5/11 


5.47 ± 0.56 


P 


0.523 ± 0.063 


0.1027 ± 0.0050 


5.5/11 


00 + 01 
u - uu -0.00 


Bf 


0.517 ± 0.044 


0.1160 ± 0.0029 


3.5/11 


13.75 ± 0.88 


B w 


0.476 ± 0.100 


0.1134 ± 0.0042 


4.1/11 


00 +0 - 05 

u - uu -0.00 


C 


0.537 ± 0.044 


0.1125 ± 0.0038 


6.3/11 


11.58 ± 0.88 



Table 3: Determination of a and a s (m z ) from fits to the first moments of the event shape 
distributions together with x 2 /d.o.f. from those fits. Also shown is the A 2 parameter from fits 
to the second moments. 
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(GeV) 


Second moments of 


1 - T 


P 


Bx 


Bw 


c 


30-50 


.0143 ± .0009 ± .0015 


.0080 ± .0006 ± .0005 


.0236 ± .0009 ± .0005 


.0104 ± .0005 ± .0005 


.1726 ± .0078 ± .0115 


50-60 


.0109 ± .0008 ± .0006 


.0063 ± .0005 ± .0008 


.0187 ± .0008 ± .0012 


.0086 ± .0005 ± .0006 


.1308 ± .0066 ± .0063 


60-70 


.0109 ± .0006 ± .0010 


.0060 ± .0004 ± .0011 


.0187 ± .0006 ± .0022 


.0088 ± .0003 ± .0013 


.1308 ± .0050 ± .0164 


70-80 


.0093 ± .0004 ± .0010 


.0053 ± .0002 ± .0007 


.0172 ± .0005 ± .0014 


.0081 ± .0003 ± .0008 


.1176 ± .0037 ± .0117 


80-84 


.0086 ± .0005 ± .0010 


.0052 ± .0003 ± .0007 


.0160 ± .0006 ± .0015 


.0081 ± .0003 ± .0008 


.1110 ± .0047 ± .0125 


84-86 


.0086 ± .0006 ± .0020 


.0054 ± .0004 ± .0014 


.0158 ± .0007 ± .0022 


.0082 ± .0004 ± .0018 


.1115 ± .0058 ± .0195 


91.2 


.0077 ± .0001 ± .0003 


.0053 ± .0001 ± .0002 


.0158 ± .0001 ± .0003 


.0076 ± .0001 ± .0002 


.1034 ± .0003 ± .0031 


130 


.0064 ± .0005 ± .0002 


.0041 ± .0003 ± .0001 


.0131 ± .0006 ± .0002 


.0069 ± .0004 ± .0001 


.0848 ± .0050 ± .0025 


136 


.0080 ± .0008 ± .0007 


.0045 ± .0004 ± .0001 


.0141 ± .0008 ± .0004 


.0076 ± .0005 ± .0002 


.0938 ± .0064 ± .0017 


161 


.0059 ± .0007 ± .0002 


.0040 ± .0004 ± .0001 


.0121 ± .0008 ± .0004 


.0070 ± .0005 ± .0002 


.0757 ± .0064 ± .0019 


172 


.0064 ± .0009 ± .0005 


.0040 ± .0005 ± .0003 


.0136 ± .0014 ± .0013 


.0068 ± .0006 ± .0005 


.0979 ± .0133 ± .0129 


183 


.0064 ± .0005 ± .0001 


.0042 ± .0003 ± .0002 


.0121 ± .0006 ± .0003 


.0067 ± .0003 ± .0002 


.0804 ± .0051 ± .0032 


189 


.0064 ± .0004 ± .0004 


.0043 ± .0002 ± .0002 


.0121 ± .0004 ± .0005 


.0071 ± .0002 ± .0002 


.0794 ± .0032 ± .0038 



Table 4: Second moments of the five event shape variables at different energy points. The two errors are respectively statistical and 
systematic. 





(1 — T) 



r 


1)1 




C 


Fit Range 


0.00-0.30 


0.00-0.20 


0.00-0.25 


0.00-0.20 


0.05-0.50 


rv (~\ 30 PpV^ 
ctg^lOU \j ( V J 


n 11 3d 
u. i ioy 


n 1 1 


nut;') 


n i nfi^ 

U. lUOO 


U.1101 


a t" l t;t 1 Pfl 1 PVTOT" 
kJ LCI L lo LILjdl KjL L Ul 


+0 0036 


±0 0034 


±0 0027 

— I—VJ . KJKJ^l I 


±0 0097 

—I— KJ . kjkJZj 1 


±0 0036 

—1—kJ. kJkJkjkj 


Systematic error 


±0.0028 


±0.0029 


±0.0016 


±0.0015 


±0.0018 


Ovprall pvnprrmp'ntal prroT" 

V/ V V_/l Oil CA Uvl llllC'llUUil C/l L Ul 


±0 0046 


±0 0045 


±0.0031 


±0 0031 

_i— kj • yjyjtj i 


±0 0040 

—1— KJ • KJKJ^KJ 


Overall theoretical error 


±0.0056 


±0.0038 


±0.0062 


±0.0088 


±0.0066 


y 2 /dof 


6.9 / 10 


8.4 / 9 


9.1 / 11 


12.0 / 12 


8.5 / 8 


Ctg^lOU ucv J 




n 1 1 1 9 


n 1 1 41 


n i da^ 


n 1 nso 


Statistical prrnr 
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—i—\J lUUU 1 
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Systematic error 


±0.0024 


±0.0013 


±0.0010 


±0.0026 


±0.0020 


Ovprall pvnpriinpntfll prror 
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Overall theoretical error 


±0.0060 
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±0.0078 
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Y 2 /d.o.f. 


10.2 / 9 


11.4 / 13 
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U. lUlo 


n 1 m 9 


n 1 1 m 
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n 1 n^9 


U. lU4io 
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Systematic error 


±0.0022 


±0.0022 
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±0.0044 
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V^/ V .L Clill W\. l_/V_^l lllll^llUCll v>l 1 KJ1. 
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Overall theoretical error 
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±0.0057 
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5.7 / 13 
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ct s yl i L vi e V J 
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±0 0026 

— 1 — KJ • UUiJ KJ 


±0 0025 

— 1 — VI • KJ KJ i^*J 


±0 0021 


±0 0016 

— 1 — KJ •KJKJ ±KJ 
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Overall theoretical error 
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6.4 / 13 
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U. 1 1UO 


n 1 1 ns 


n 1 1 1 a 


n i 

U. 1UOO 


n 1 1 1 s 

U.lllO 


O IdLlo LlUdl tilUl 


+D DDI 4 


XU.UUIO 


+n nm 1 

xu.uui ± 




+n nm 8 


Systematic error 


±0.0012 


±0.0010 


±0.0014 


±0.0012 


±0.0014 


Overall experimental error 


±0.0018 


±0.0016 


±0.0018 


±0.0016 


±0.0023 


Overall theoretical error 


±0.0057 


±0.0033 


±0.0067 


±0.0078 


±0.0055 


xVcl.o.f. 


4.4 / 11 


8.2 / 13 


28.0 / 12 


10.6 / 13 


5.7 / 8 



Table 5: a s measured at ^/s = 130, 136, 161, 172, 183 and 189 GeV from fits of the event shape 
variables to theoretical predictions with combined fixed order and resummed calculations. The 
fit ranges, the estimated experimental and theoretical errors and the fit quality are also given. 
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(GeV) 


a s (from 


T, p, B T 


, B w ) 




(from T, p, B T , B w , C) 


30-50 





1400 


± 


0.0056 


± 


0.0107 






50-60 





1260 


± 


0.0073 


± 


0.0088 






60-70 





1340 


± 


0.0060 


± 


0.0087 






70-80 





1210 


± 


0.0064 


± 


0.0082 






80-84 





1200 


± 


0.0057 


± 


0.0089 






84-86 





1160 


± 


0.0061 


± 


0.0082 






91.2 





1221 


± 


0.0020 


± 


0.0066 






130 





1122 


± 


0.0038 


± 


0.0060 





1128 ± 0.0038 ± 0.0063 


136 





1116 


± 


0.0042 


± 


0.0060 





1111 ± 0.0043 ± 0.0061 


161 





1041 


± 


0.0052 


± 


0.0054 





1041 ± 0.0054 ± 0.0054 


172 





1075 


± 


0.0054 


± 


0.0056 





1084 ± 0.0056 ± 0.0055 


183 





1089 


± 


0.0022 


± 


0.0056 





1088 ± 0.0023 ± 0.0055 


189 





1106 


± 


0.0017 


± 


0.0058 





1105 ± 0.0018 ± 0.0058 



Table 6: Summary of a s values as determined from event shape variables at different centre-of- 
mass energies. The a s values for -^/i < rriz were determined JQE] only from four event shape 
variables for which analytical calculations were available at that time. 
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Figure 1: Distributions for thrust, T, scaled heavy jet mass, p, total and wide jet broadenings, 
Bt and B\y, and the C-parameter at y/s = 189 GeV in comparison with QCD model predictions. 
The errors shown are statistical only. 
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Figure 2: The first moments of the five event shape variables, 1 — T, p, B T , B w and C, as a 
function of the centre-of-mass energy, compared with several QCD models. 
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Figure 3: The first moments of the five event shape variables, 1 — T, p, B?, Byy and C compared 
to the results of a fit including perturbative and power law contributions. 
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Figure 4: The second moments of the five event shape variables, 1 — T, p, B?, B\y and C 
compared to the results of a fit including perturbative and power law contributions. The 
parameters a and a s are fixed to the values obtained by the fits to the first moments. The 
A 2 /s term is negligibly small for p and Bw but is necessary to reproduce the behaviour of 1 — T, 
B T and C. 
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Figure 5: a) a s measurements from event shape distributions as a function of the centre-of-mass 
energy. The errors shown are experimental only. The solid and dashed lines are fits with the 
energy dependence of a s as expected from QCD and with constant a s , respectively, 
b) a s values as determined by L3 from the r lifetime measurement, Z lineshape and event shape 
distributions. The line is a fit to the QCD evolution function to the measurements made from 
event shape variables. 
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